A simple method for the determination of slowly varying refractive index profiles 
from in situ spectrophotometric measurements 
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Reliable control of the deposition process of optical films and coatings frequently requires mon- 
itoring of the refractive index profile throughout the layer. In the present work a simple in situ 
approach is proposed which uses a WKBJ matrix representation of the optical transfer function of a 
single thin film on a substrate. Mathematical expressions are developed which represent the minima 
and maxima envelopes of the curves transmittance-vs-time and reflectance-vs-time. The refractive 
index and extinction coefficient depth profiles of different films are calculated from simulated spectra 
as well as from experimental data obtained during PECVD of silicon-compound films. Variation of 
the deposition rate with time is also evaluated from the position of the spectra extrema as a function 
of time. The physical and mathematical limitations of the method are discussed. 
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I. INTRODUCTION 

In optical coatings, two different kinds of graded index 
profiles can occur: (i) intentionally graded profile, when 
precisely controlled inhomogeneous optical coatings (for 
example, rugate filters) are produced;!!! (ii) accidentally 
created inhomogeneities that may arise from the insta- 
bility of the deposition process. The latter may be due 
to different factors ,00 including the change of microstruc- 
ture of the film deposited, the contamination, the lack of 
control of the fabrication parameters, or the interaction 
of the growing film with the substrate. 

In order to control the film quality it is, therefore, 
very important to calculate the refractive index depth 
profiles from in situ measurements. Bovard,Q inspired by 
the popular transmittance envelope methocLfox the calcu- 
lation of refractive index dispersion curvesJjLj developed 
an approximative method based on the transmittance en- 
velopes of the curve transmittance-vs-time. In this work, 
we describe a generalized reflectance and transmittance 
envelope method based on WKBJ approximation. A ma- 
trix representation of thin films with slowly varying index 
of refraction is used to obtain maxima and minima enve- 
lope expressions, which are then analytically solved and 
the refractive index profile is extracted. We apply this 
method to evaluate the uniformity of the refractive index 
in amorphous hydrogenated silicon nitride (SiNi. 3 ) films 
deposited in low pressure plasma. 



II. EXPERIMENTAL APPARATUS 

Before we present the theoretical background of our 
calculations, we first describe the experimental method- 



ology for the preparation of optical films. The films stud- 
ied in this work were grown on a radiofrequency (13.56 
MHz) -powered electrode (18 cmin diameter) in a plasma 
system (Fig. |l|) described earlier H A mixture of silane and 
ammonia (typically 1:3 ratio) was used for the deposition 
of SiNi.3 on both glass and silicon wafer substrates, us- 
ing a working pressure of 40 mTorr. In order to ease the 
optical reflectance analysis, the back surface of the glass 
substrate was roughened. 

During the optical reflectance experiment, a white light 
beam (100 W halogen lamp, Oriel) hit the surface of 
the growing film at an angle of incidence close to the 
normal. The reflected light was collected by an opti- 
cal fiber and brought to a spectrometer (Multispec I, 
Oriel) equipped with a photodiode array (1024 diodes, 
Instaspec II, Oriel). The acquisition time of the detector 
was kept between 0.3 s (Si substrate) and 0.6 s (glass 
substrate). Light from the plasma was recorded as a 
background and subtracted from the spectra. Each 4-6 
seconds during the deposition, a spectrum was recorded. 
Example of reflectance data measured in situ with the 
optical monitor is shown in Fig. |^. One can see the vari- 
ation of the maxima and minima of reflectance both with 
time and wavelength. A cross section of the spectrum 
in Fig. |^(a) is shown in Fig. [|(b), which represents the 
variation with time of the reflectance at a wavelength of 
789 nm. 

In addition, ex situ reflectance and spectroellipsomet- 
ric measurements have been made using a Perkin-Elmer 
lambda- 19 spectrophotometer and a variable angle spec- 
troscopic ellipsometer (VASE, J. A. Woollam Co.), re- 
spectively. 
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III. DESCRIPTION OF THE METHOD 



A. Refractive index profile 



The method used to determine the refractive index of 
the film at any time of the deposition process is based on 
the evaluation of the envelope curves connecting the re- 
flectance minima and maxima (dashed lines in Fig. 0(b)). 
Following a WKBJ (Wentzel-Kramers-Brillouin-Jeffries) 
approximation^ one can express the "optical transfer 
function" of an inhomogeneous thin film by the following 
characteristic matrix£3 



M = 



■ cos 5 



■ sin S 



V"m n out sin 5 



T^out 



cos 5 



(1) 



with S = (2ir/\) J [n{z) — \k{z)\dz for normal incidence, 
where ni n and n out are the film refractive index values 
near the film-substrate interface and near the surface, re- 
spectively, d is the thickness of the layer, z is the distance 
from the interface inside the layer, and n(z) is the refrac- 
tive index profile in the layer. From this characteristic 
matrix one can obtain expressions for the transmittancc 
T and reflectance R of a film on a substrate: 



n n in 



Untout exp(<5) 



R = 



rout hn exp(2<5) 



cxp(25) 



1 + r out T in exp(2(S) 



(2) 



(3) 



Here, t out = 2n /(n out + n ), t in = 2n in /{m 
rout = (n - n out )/(n Q + n out ) and r in = (n in - n s + 
ik s )/(rii n + n s — ik s ). The sign over Ti n indicates that 
this parameter can have both complex or real values, as 
one can consider the reflectance problem for both trans- 
parent and absorbing substrates. Not surprisingly, the 
expressions above are very similar to those describing the 
homogeneous case, and they translate very well the fact 
that the WKBJ approximation is equivalent to taking 
into account only the reflections occurring at the inter- 
faces in the coatings, while the particular shape of the 
refractive index inside the layer is neglected. 

From the approximate expressions (Q) and (^) , one can 
obtain expressions for the minima and maxima of the 
transmittance or reflectance curves, T(z) or R(z). From 
the data available (monitored R(z) or T(z) values) one 
can then evaluate experimentally the envelopes and solve 
the envelopes expressions (such as Eqs. @j and (g) , see 
below) for n(z). However, several assumptions must be 
respected: (i) the value of the film extinction coefficient is 
low (k <C n) , so that only the k appearing in the exponen- 
tial terms is significant; (ii) once a material is deposited, 
its refractive index value is not altered by the following 
deposition on top of it; (iii) the substrate is considered 
semi-infinite, which means that the contribution to the 
reflectance from the back side is negligible. 



1. Transmittance 



BovarcH has developed a method for the transmit- 
tance which consists of solving the equations for the 
envelopes of the transmittance minima and maxima 
for n jn , n ou j and the absorbance parameter A = 

exp[(— 4ir/\) J Q d k(z)dz\. In the following we show a gen- 
eralized solution, for low or high index substrates, based 
on his previous workJa For d = 0, one obtains: 



N t + r n /N?-n 2 n 2 s 



1/2 



Ni 



2enon s 



(4) 



and for d > 0: 



2n in n s n e (T, 



T 

± rr 



n 



(ni n - n 2 a Y 



T —T 

1 max 1 mm 



T, 



max- 1 - nun 



A(z) = 



£ y^max/^-rr 



T infant \fT max jT m i n + 1 



1/2 

(5) 
(6) 

and 



with e = —1 for n,;„ > |n s |, e = +1 for rii n < \h s 
n = 1 for nj n > n \h s \, r] = -1 for nf n < n \n s \. 



2. Reflectance 



The solution to the inverse problem, i.e. finding the 
variation of the optical parameters with depth from 
Eq. ([?]), consists of solving the following minima and 
maxima envelopes relations for three unknowns Ti n , v out 
and A: 



Rn 



r 2 

out 



2er ou t$t{r in )A+ \f in \ 2 A 2 



1 - 2er ou M{r m )A + r 2 out \r in \ 2 A 2 ' 



Rn'i.ax ~ 



2er n 



2 A 2 



^ 'out I ^*in 

\ 2 A 2 ' 



(7) 
(8) 



Here, 5?(rj n ) is the real part of r in . The system of Eqs. (^) 
and (||) offers only two equations. Therefore, to solve 
it we assume that at the first stage of the deposition 
(d — * 0), the film is homogeneous (n out — * xii n ) and trans- 
parent (A — > 1). Using this approximation, and solving 
for rLj n , we find 



2 



where 



_ (e + 1)-Rmin — ( e — 1)-Rmax + 2 

(e - l)i? maa ; — (e + l)R mm + 2 

(/ n 2 j^i \ Keeping the value of n.j„ constant, and solving for r ou t 

N r + m / N? 2 , (9) and for A when d > 0, the calculated solutions are as 

V Hs J follows: 



r out{ z ) — ' 

A(z) = 1/2 



|(|-CC-if4 + (f- W 



1/2 



( ^ ^o-ui ) {Rmax Rmin } 



C*oui^(^m) [^otit i^RminR' 



max ■ lL max 



Rn 



mzn ±L max 



Rn 



(10) 

(11) 



where 



|^*m| [Rmax Rmin} ? 
5 23^(7*^) (2 Rmax Rmin} (2 Rmin Rmax Rmin Rmax}-} 



B 2 D 



D 2|?"m| [Rmax Rmin} H~~ 8^(^*?'n) [(-^mas H~~ Rmin) ~1~ 2(1 Rmax Rmin} 



-\- < lR rrlau: R, 



max l 'min 



1 1 



Rmin R' 



- 1 



mm ll max 



and 



c= < 



> 



-1 if 



< 



B 

T 



The refractive index profile n out (z) is then calculated 
from n out = n (l - r out )/(l + r out ). 

To illustrate the method with an example, we ap- 
plied the above expressions to the envelopes surround- 
ing the simulated spectra of a l-/^m-thick transparent 
(k = 0) layer with a linear profile and a deposition rate of 
60 nm/min (Fig. ||). The profiles are expressed as a func- 
tion of time instead of distance z, to simulate the type 
of data obtained experimentally with an optical moni- 
tor. The calculated profile in Fig. |3|(b), although clearly 
approximate, is in very good agreement with the model 
linear profile used in the simulation. 



B. Extinction coefficient profile 

Calculation of k(z) from *4(z) in Eq. ([j|) or ( [Tl| ) is done 
using the following relation: 



k(z) 



-X d 

47T dz 



\nA{z) 



(12) 



Figure || shows the ^4(z) profile obtained from the en- 
velopes of Fig. ||(a). The corresponding extinction coeffi- 
cient profile is shown in Fig. |](b). One can see that small 
deviations of A(z) from 1 (the value it should have) have 
strong impact on the value of the calculated extinction 
coefficient, which may reach unexpected negative values. 
Therefore, experimental precision of the envelopes must 
be very high if one is to find reliable k values. 



C. Deposition rate and physical position in the layer 

It is interesting to express the results in Figs. |^ and || 
as a function of physical position inside the layer instead 
of deposition time. When the deposition rate is constant 
during the deposition, the transformation of the results 
as a function of thickness is obvious. Unfortunately, it 
may not be the case for numerous experiments. If the 
deposition rate varies during the deposition, one can still 
use the fact that the optical thickness between two suc- 
cessive extrema is equal to a quarterwave. Using the 
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relation mA/4 = J Q n(z)dz (m is an integer representing 
the interference order) and assuming that the deposition 
rate is constant between successive extrema, one obtains 
the following relation: 

i^r-Lf f\ m , (is) 

4 *2 - tl J tl 

where zi, Z2 and tx, t2 correspond to the positions of two 
successive extrema on the thickness (distance) and time 
scales, respectively. The relationship between the time 
and distance scales can be obtained from ti = (zi = 0), 
where the first extremum occurs. Figure |^(a) shows the 
calculated evolution of the deposition rate from the sim- 
ulation in Fig. m The deposition rate varies around the 
expected value (60 nm/min) within an interval of plus 
or minus 4 nm/min. One can hence calculate the rela- 
tion between the deposition time axis and the distance in 



the layer (see Fig. |](b)). It appears that variation of the 
calculated deposition rate around its overall mean value 
is negligible and that total thickness found is consistent 
with the expected value of 1 /xm. 



D. Error calculus 

1 . Transmittance 

A good estimate of the error in determining the n(z) 
and A(z) values may be obtained from the calculation of 
the derivatives of Eqs. (^)-(^|). In the case of the trans- 
mittance, assuming that n s is known precisely and that 
the error of T max and T m „ are the same and not related, 
we have 



2(n n s 



T 



max \ 1 mm 

T — T / T 

± m ax - 1 - mm / 



AT 



(14) 



n ut 



n C„ 



with 



Similarly: 



T 



"7 



An*. 



A*A. ^ i ^\/T rriax T rn i n AT ( 2,nQn OU f An ou t 

•A Tmax Train T 



n out 



7ig T^out 



(15) 



(16) 



(17) 



2. Reflectance 

Expressions for the error in the envelope reflectance 
method can be found from the derivation of Eqs. 
( pd| ) as functions of R and n in . The resulting equations 
are not shown here, but have been used for the determi- 
nation of the error bars in the figures appearing in this 
work (see Figs. |, [|, 0, [| and ||). The error on R max 
and R m i„ have been fixed to a value of AR =1% and 
corresponds to the experimental error on the measured 
extrema. This value doesn't take into account the error 
generated during the calculation of the envelopes. 



IV. APPLICATION TO EXPERIMENTAL DATA 

The method has been tested experimentally for silicon 
nitride films deposited on both transparent (glass) and 
absorbing (silicon) substrates. The results are shown, 
respectively, in Figs. ^ and 0. In these experiments, the 



deposition conditions were deliberately not precisely con- 
trolled: as a result, a refractive index gradient in the 
layers is observed. The profile in Fig. || is realistic; The 
variation in n(z) appearing in Fig. |^ may be mostly due 
to oxygen contamination by water vapor present in the 
earlier stage of the film growth. Additionally, n(z) in 
Fig. H may vary because of unstable deposition parame- 
ters, such as oxygen leak. Index variations could also be 
due to changes in the microstructure of the layer (varying 
porosity). Figures ^|(c) and|7](c) show a higher deposition 
rate value in the initial stage of the deposition process; 
further investigations are on the way to identify further 
experimental evidence of these phenomena. 

Repeating the calculation for several wavelengths, the 
value of the refractive index dispersion at different po- 
sitions in the layer can be obtained, as shown in Fig. ||. 
The distance has been calculated from Eq. (p"3|); the mean 
value of the deposition rate has been used (22 nm/min). 
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V. DISCUSSION 

A. Experimental error 

One must be very careful when using an optical mon- 
itor for the quasi-continuous measurement of transmit- 
tance or reflectance with time in a plasma environment. 
For measurements at a single wavelength in a glow dis- 
charge, use of a lock-in detection system would help to 
subtract the light generated by plasma. In our case, with 
a photodiode array, such light chopping was not used. In 
addition, it was found that the main source of experimen- 
tal error was the reflection and the scattering of light on 
different metallic components in the deposition chamber. 
Those were suppressed by installing a diaphragm between 
the optical monitor and the sample, and by blackening 
the surface of the sample holder. Such perturbations will 
hardly change the phase information of the signal, but 
they can have a significant effect on the signal amplitude. 



B. WKBJ approximation limitations 



As pointed out in Sec. Ill, the matrix representation 
in Eq. ([!]) follows the WKBJ approximation, which holds 
when 



n 2 (z) 



< 



(18) 



The method has been tested on simulated envelopes ob- 
tained from a slightly more complex refractive index pro- 
file. Figure ||(a) shows the model profile used and the cal- 
culated profile, while Fig. p|(b) shows an evaluation of the 
left (solid line) and right (dashed line) part of Eq. (18). 
One can see that the relation Jl8| ) holds throughout the 
profile. Consequently, the WKBJ approximation is valid 
for these profiles, like it is for the other less complex lin- 
ear profiles studied in this work. 



C. Envelopes-related limitations 

Precise determination of the envelopes from experi- 
mental spectra is delicate. Algorithms for finding en- 
velopes were developed and are used for the determi- 
nation of refractive index dispersion from ex situ spec- 
trophotometric measurements. OE3 In these cases, the 
spectra and its envelopes are functions of wavelength in- 
stead of time or distance, so that limit conditions at the 
right and left sides of the spectra can be used in order 
to improve the accuracy of the envelope computed. No 
such limit conditions exist in our case. Consequently, 
simple cubic or linear spline has been chosen to compute 
the envelopes. As a result, the left and right extremes 
of the spectra are more likely to generate errors in the 
calculus. In addition, to increase the accuracy of the en- 
velopes, the layer must reach a minimum thickness, so 



that the spectra show at least 3 extrema. The value of 
this minimum thickness will depend on the index of the 
film and substrate materials used, and on the wavelength 
of the probing light. The shorter the wavelength, the 
smaller the minimum thickness will be. Once such min- 
imum thickness is reached, the method can be used to 
follow, in real time, the evolution of the refractive index 
profile during the film growth. 



D. Change of sign of £ within the layer 



In Sec. Ill, the sign parameters e, ?y and £ can all reach 
values of +1 or -1, depending on the different values of 
the refractive index of the substrate and of the film. The 
former two are easy to determine, as they depend only on 
the values of n s and Ui n , and do not depend on z. The 
latter one (() can be more difficult to determine since 
its value depends on n OMt (z) and is changing along the 
profile. Figure 10 illustrates the effect of a change in 
the sign of C on the refractive index profile: Similarly 
to Fig. ||, a linear profile (dashed line) has been used to 
simulate the variation with time of the reflectance during 
the deposition of the graded layer. The envelope method 
developed above has been used to calculate the refrac- 
tive index profile twice, i.e. with £ value of +1 and -1, 
respectively. Around t « 7 min., the variation of r out , 
Rmin and Umax is such that the value of C in Eq. jltj| ) 
should be changed from -1 to +1 in order to calculate an 
accurate refractive index profile. 



E. Ex situ analysis 

For comparison additional ex situ spectrophotomet- 
ric and spectroellipsometric analysis of samples from 
Sec. ^ have been performed. Using the n(z) profile 
for SiNi.3 layer on glass calculated with the envelope 
method (Fig. ^|), cutting this profile in small sublayers, 
and introducing a small n(A) dispersion, it was found 
that the reflectance calculated from this profile repro- 
duce well the ex situ measured reflectance (Fig. |ll|(b)) 
and the in situ measured reflectance (Fig. 0(c)). In ad- 
dition, a simple model considering a single homogeneous 
layer (dotted line in Fig. [n](a)) has been used to repro- 
duce more precisely the ex situ reflectance (Fig. |ll](b)); it 
was found that this model gives almost the same in situ 
reflectance variation with time (Fig. 11(c)) as the graded 
profile (solid line in Fig. |ll](a)). However, the ex situ 
spectra provide very little information about the exact 
index profile; in fact, different changes in the thickness 
and the n(z) values of the graded profile in Fig. |Tl](a) 
could also result in a good match with the ex situ data 
in Fig. |TT](b) . Therefore, ex situ measurements do not 
offer enough reliable informations about the real profile 
n(z). Ex situ spectroellipsometric measurements led to 
the same observations. One can conclude that: (i) cx 
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situ analysis can give a false impression of homogeneity, 
due to the fact that several different models can fit ex 
situ data (there are several local minima of the merit 
function used in the fitting procedure); (ii) in situ anal- 
ysis using the envelope method proposed in this work 
avoids fitting problems, but small errors in the exper- 
imental data can lead to large uncertainties in the n(z) 
and k(z) profiles. The advantage of in situ measurements 
is clear ("4i.'s easier to investigate a murder if you can see 
it in situ"Ej), but finding n(z) very precisely requires very 
sensitive measurements. 



VI. CONCLUSION 

An envelope method has been developed for the de- 
termination of refractive index profile n(z) from in situ 
spectrophotometric measurements (reflectance or trans- 
mittance). The method is based on the WKBJ approxi- 
mation valid for slowly varying refractive index profiles. 
Generalized analytical expressions, suitable for high or 
low index substrates, are given for the calculation of n(z) 
from in situ transmittance measurements. A more gener- 
alized solution, valid also for non-transparent substrates 
was calculated and it is used for the determination of 
n(z) from in situ reflectance measurements. The instan- 
taneous deposition rate and its variation during the depo- 
sition process are derived from the calculated n(z) profile. 
The main feature of the envelope method developed in 
this work is its analytical character: Unlike usual in situ 
spectrophotometric and ellipsometric analysis, no non- 
linear fit algorithms or "initial guess" of the refractive 
index has to be used in order to ascertain n(z). In fact, 
the simplicity of this envelope method makes it attractive 
for rapid first approximation analysis to generate initial 
data for more sophisticated methods such as in situ el- 
lipsometry. 
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FIG. 2. (a) Typical spectra evolution measured with the 
optical monitor, (b) Closer look at the reflectance modula- 
tion with time at A = 789 nm. 
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FIG. 4. Linear refractive index profile: (a) A(t) profile, 
and (b) k(t) profile (dashed line: model, solid line: calcu- 
lated) . 



FIG. 5. Linear refractive index profile: (a) deposition 
rate values between successive extrema, (b) physical posi- 
tion in the layer as a function of deposition time (dashed 
line: model, solid line: calculated). 
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FIG. 7. Experimental data for an inhomogeneous SiNi.3 
film on silicon: (a) Reflectance evolution, (b) calculated n(z) 
profile, (c) calculated deposition rate variation. (A = 826 nm) 




FIG. 8. Experimental data for an inhomogeneous SiNi.3 
film on glass: Evolution of the refractive index dispersion 
within the layer (same film as in Fig. ^|). 
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FIG. 9. Half-period-sinus refractive index profile: (a) 
model (dashed line) and calculated (solid line) n(z) profiles, 
(b) WKBJ validity condition from Eq. ( |l8| ) (dashed line: 
1/A, A = 800 nm). 
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FIG. 11. Ex situ optical analysis of S1N1.3 on glass (same 
sample as in Fig. ^): (a) n(z) profile calculated using the 
envelope method (solid line, same as 

Fig- ib)) 

and ho- 
mogeneous one-layer model (dotted line) used to fit the 
ex situ measurements (in (b)); (b) ex situ reflectance data 
(gray line), and reflectance calculated using the profiles in 
(a) (black solid line: graded profile; dotted line: optimized 
homogeneous profile); (c) in situ reflectance measurement 
(gray line), and reflectance calculated using the profiles in 
(a) (black solid line: graded profile; dotted line: optimized 
homogeneous profile). 
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FIG. 10. Effect of of £ on n(t) profiles (dashed line: model 
profile, solid: calculated profiles, A = 800 nm). 
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